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RESEARCH  FINDINGS: 

This  research  uses  a  simple 
experimental  chemical  model  to  understand 
the  nonlinear  dynamics  of  the  complex 
biochemical  network  of  hemostasis  (blood 
clotting).  The  chemical  model  is  based  on  a 
modular  mechanism  describing  initiation  of 
clotting.  Briefly,  we  constructed  a  chemical 
model  of  the  hemostasis  network,  which 
contains  approximately  80  reactions.  Our 
model  contained  three  interacting  modules; 
each  module  consisted  of  a  single  chemical 
reaction  with  appropriate  kinetics.  In  the 
model,  a  threshold  response  was  regulated 
by  the  interaction  of  the  three  modules. 
Previously  published  work  describes  this 
approach  in  detail  ( PNAS  2006  103:15747- 
15752).  Using  this  model  system,  we 
demonstrated  that  the  spatial  localization  of 
bacteria  profoundly  affected  coagulation  of 
human  and  mouse  blood  and  plasma. 

Spatial  Localization  of  B.  cereus  controls 
coagulation 

We  hypothesized  that  initiation  of 
coagulation  by  bacteria  would  be  regulated 
by  the  spatial  localization,  not  the  total 
amount,  of  bacteria  (Fig  la).  In  other 
words,  for  bacteria  that  activate  coagulation 
factors,  coagulation  would  only  occur  when 
a  cluster  of  bacteria  forms.  This  hypothesis 
was  based  on  previous  experiments  with 
human  blood  and  plasma  that  showed  i) 
stimuli  must  exceed  a  local  threshold 

.17 

concentration  to  initiate  coagulation,  *  and 
ii)  this  threshold  response  to  concentration 
leads  to  a  spatial  threshold  response,  in 
which  coagulation  initiated  on  a  patch  of 
stimulus  above,  but  not  below,  a  threshold 


Threshold  for  initiation  of  coagulation 
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Figure  1 .  Human  blood  plasma  coagulates  on 
spatially  localized  B.  cereus  (a)  A  schematic 
drawing  illustrates  the  coagulation  potential  of  a 
single  bacterium,  a  single  small  surface  patch  of 
bacteria,  surface  patches  of  bacteria  that  are 
separated,  and  surface  patches  of  bacteria  that  are 
clustered  together,  (b)  A  chart  compares  the 
clotting  times  of  human  blood  plasma  on  the  same 
amount  of  bacteria  either  clustered  in  a  large  patch 
or  dispersed  in  solution,  (c  and  d)  Fluorescent 
photographs  of  a  microfluidic  chamber  used  to  test 
coagulation  (blue  fluorescence)  of  human  blood 
plasma  on  clustered  bacteria  (green  fluorescence) 
patterned  as  patches  of  different  sizes. 

Coagulation  occurs  rapidly  on  a  large  patch  (d),  but 
not  on  an  array  of  smaller  patches  containing  the 
same  amount  of  bacteria  (c). 


To  test  this  hypothesis,  we  compared  the  clot  time  of  human  blood  plasma  exposed  to 
bacteria  dispersed  in  solution  to  bacteria  clustered  on  the  surface  of  a  microfluidic  chamber 
Bacillus  cereus  (B.  cereus )  spatially  localized  to  a  surface  cluster  rapidly  initiated  coagulation. 
However.  B.  cereus  dispersed  in  solution  at  concentrations  of  up  to  107  colony-forming  units 
(CFU)  mL'1  (Fig  lb).  In  a  second  experiment,  we  used  microfluidics5’6  and  micropattemed 
surfaces  to  control  the  spatial  distribution  of  bacteria  and  to  demonstrate  that  the  size  of  the 


cluster,  rather  than  amount  of  bacteria,  can  control  the  rate  of  initiation  of  coagulation  of  human 
blood  plasma.  On  smaller  patches  (90  Dm)  spaced  far  apart  (400  Dm),  coagulation  was  slow, 
initiating  on  the  first  patch  in  9  min  □  1  with  clotting  on  all  the  patches  in  the  array  in  22  □  3 
min  (mean  □  S.E.),  indicating  that  the  individual,  90  Dm  patches  were  below  the  size  necessary 
to  initiate  coagulation  rapidly  (Fig.  lc).  However,  when  the  same  number  of  bacteria  were 
patterned  closer  together  to  form  a  large  patch,  coagulation  initiated  rapidly  in  5  □  1  min  (mean 
S.E.)  over  the  entire  patch  (Fig.  Id)  (p-value  <  0.01  in  comparison  to  initiation  on  the  first  90 
□m  patch,  and  <  0.005  in  comparison  to  initiation  on  the  entire  set  of  patches). 

B.  cereus  initiates  coagulation  of  flow  ing  whole  blood 

To  test  if  B.  cereus  would 
initiate  coagulation  in  the 
presence  of  flow,  human 
whole  blood  was  flowed 
over  localized  colonies  of 
B.  cereus  in  microfluidic 
channels  (Fig.  2).  Clusters 
of  bacteria  in  microfluidic 
channels  were  made  by 
encapsulating  bacteria  in 
gel  microdroplets  (GMDs).7 
GMDs  consisted  of 
colonies  of  bacteria  and 
magnetic  particles  1  pm  in 
diameter  contained  in 
agarose  spheres 
approximately  50  pm  in 
diameter;  the  magnetic 
particles  allowed  the  GMDs 
to  be  trapped  in  the 
microfluidic  channels  by  a 
magnet  incorporated  in  the 
device  near  the  channel 
(Fig.  2a).8  Clusters  of  B. 
cereus  initiated  coagulation 
of  flowing  human  whole 
blood  in  3-13  min  (Fig.  2b,c),  whereas  coagulation  did  not  occur  until  48-59  min  in  experiments 
with  the  control  strain  of  E.  coli  (Fig.  2d.  p  <  0.001). 


Figure  2.  Human  whole  blood  coagulated  on  spatially  localized  B. 
cereus  in  the  presence  of  flow,  (a)  A  microphotograph  shows  the 
microfluidic  device  used  to  flow  human  whole  blood  over  colonies  of 
bacteria.  The  bacteria  were  localized  in  gel  microdroplets  (GMDs)  that 
also  contained  magnetic  particles.  Magnets  were  used  to  localize  the 
GMDs  in  the  device,  (b)  A  brightfield  image  shows  coagulation  of 
human  whole  blood  on  GMDs  containing  B.  cereus  expressing  GFP 
(overlaid  green  fluorescence),  (c-d)  Microphotographs  of  fluorescence 
show  coagulation  of  human  whole  blood  (blue)  on  GMDs  containing  B. 
cereus  (green)  (c),  but  on  not  E.  coli  (green)  (d).  White  dashed  lines 
outline  the  channel  walls.  Images  were  taken  1 1  min  after  blood  was 
introduced  into  the  device,  (e)  A  graph  shows  the  clot  times  of  flowing 
whole  blood  on  colonies  of  B.  cereus  and  E.  coli  in  microfludic  devices. 
Each  data  point  represents  an  independent  experiment. 


Clusters  of  B.  anthracis  initiate  coagulation  in  mice 

We  used  a  mouse  model  to  test  whether  the  correlation  between  clustering  of  B.  anthracis  and 
local  coagulation  was  also  observed  in  vivo.  Mice  were  injected  with  B.  anthracis  vegetative 
bacteria,  and  then  lung,  heart,  spleen,  and  liver  tissues  were  harvested  quickly  (30  or  90  min) 
after  the  injection.  Rapid  harv  esting  was  used  to  minimize  initiation  of  coagulation  by  the 
immune  response.0,10  After  harvesting,  histological  sections  of  the  mouse  tissues  were  scored  for 
the  percentage  of  vessels  showing  fibrin  clots.  In  the  control  experiment,  two  mice  were  injected 
with  104  bacteria  per  mouse;  we  hypothesized  that  clusters  of  bacteria  would  be  less  likely  to 
form  at  such  a  low  dose.  In  these  control  mice  (Fig.  3a),  no  clusters  of  bacteria  and  no  fibrin 
clots  were  observed  in  any  tissue.  Then,  two  mice  were  injected  with  a  higher  dose  of  108 
bacteria  per  mouse;  we  hypothesized  that  the  formation  of  bacterial  clusters  would  be  more 
likely  at  this  high  dose.  In  the  mice  receiving  108  bacteria,  clusters  of  B.  anthracis  were 
observed  in  the  microvasculature  of  the  lungs  (Fig.  3b),  but  not  in  any  other  organs  30  min  after 
the  injection. 
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Figure  3.  Clusters  of  B.  anthracis  rapidly  initiate  coagulation  in  mice,  (a-b)  Hematoxylin  and  eosin 
(H&E)  stained  histological  sections  of  mouse  lung.  The  sections  outlined  in  black  were  obtained  with 
100*  magnification,  (a)  Pulmonary  vessels  from  control  mice  injected  with  a  low  dose  of  bacteria 
show  no  clusters  and  no  coagulation,  (b)  Pulmonary  vessels  from  mice  injected  with  a  higher  dose  of 
bacteria  show  clustering  of  bacteria.  Coagulation  in  vessels  (large  magenta  regions)  occurred  on 
clusters  of  bacteria  (chains  of  rod-like  bacteria  are  seen,  blue)  within  30  min.  Digitally  magnified 
portions  of  a  vessel  are  shown  in  images  on  the  left  and  right. 


Confinement  regulates  initiation  of  coagulation  by  tissue  factor 

Next,  we  hypothesized  that  bacteria  could  also  initiate  coagulation  if  sufficiently  limited 
mass  transfer  caused  active  factors  to  accumulate  above  the  local  threshold  concentration 
required  for  coagulation.  Confinement — enclosing  the  source  of  diffusing  molecules  within  a 
boundary  that  these  molecules  cannot  cross — is  a  simple  way  to  limit  mass  transfer.  First,  we 
tested  in  general  whether  confinement  could  affect  the  outcome  of  coagulation,  by  using 
microfluidics  to  confine  microparticles  carrying  a  classical  stimulus,  tissue  factor  (TF)  (Fig.  4). 
We  found  that  small  TF-coated  magnetic  silica  microparticles  (“TF-carrying  beads”)  did  not 
initiate  clotting  of  human  normal  pooled  plasma  for  over  40  min  when  in  a  less  confining  space 
(160  x  160  pm2  channels)  (Fig.  12?)  even  when  the  beads  were  located  close  to  a  wall.  However, 
equal-sized  or  even  slightly  smaller  TF-carrying  beads  initiated  clotting  in  under  10  min  when 
trapped  in  a  tightly  confining  space  (10  x  10  pm2  channels)  (Fig.  1C)  (N=5  beads  in  2  separate 
experiments,  p  <  0.001).  This  result  indicates  that  TF-carrying  beads  smaller  than  threshold  size 
did  not  initiate  coagulation  in  a  less  confined  environment,  but  did  initiate  coagulation  when  they 
were  highly  confined  in  a  space  with  limited  diffusion. 


Coagulation  of  human  plasma  is  a 
complex  process  that  has  been  shown  to 
exhibit  a  threshold  response  (on  and  off)  to 
the  concentration  of  soluble  active 
coagulation  factors  such  as  thrombin.11  We 
hypothesized  that  the  observed  response  to 
confinement  was  due  to  accumulation  of 
soluble  active  factors,  but  it  was  possible 
that  unknown  interactions  between  the 
plasma  and  the  device  also  played  a  role.  To 
exclude  this  possibility,  we  tested  whether 
the  response  to  confinement  could  be 
predicted  by  a  numerical  simulation  of  TF- 
initiated  blood  clotting,  using  a  set  of  40  rate 
equations  describing  enzyme  kinetics  and 
stoichiometric  inhibition  in  blood 
coagulation  without  surface  interactions12'14. 
We  used  a  simulation  with  geometry  similar 
to  that  of  the  experimental  microfluidic 
device  and  simulated  the  response  of  plasma 
to  small  regions,  or  patches,  of  1  nM  TF  in 
either  10  pm  or  160  pm  wide  channels  (Fig. 
5).  This  simulation  was  conducted  in  2- 
dimensional  space  (2D)  to  reduce 
computational  demands.  The  clot  time  was 
defined  as  the  time  at  winch  fibrin  reached 
half  of  its  maximum  concentration  (i.e.,  at 
4.15  pM  fibrin).  As  predicted,  patches  of  TF 
initiated  the  coagulation  network  quickly  (~ 
8  min)  in  confining  channels,  but  not  in  the 
less  confining  channels  (>  60  min).  This 
finding  agrees  with  the  experimental  results 
(Fig.  4)  and  confirms  that  coagulation 
kinetics  and  spatial  constraints  alone  are 
sufficient  to  achieve  the  confinement  effect. 
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Figure  4.  Initiation  of  blood  coagulation  by  TF-carrying 
beads  is  regulated  by  environmental  confinement.  (A)  A 
schematic  drawing  shows  the  geometry  of  the 
experimental  device  containing  channels  of  different 
sizes  (cross  sections  of  160  x  160  pm2,  40  x  40  pm2, 
and  10x10  pm2)  and  TF-carrying  beads  (Red).  (B,C) 
Initiation  of  coagulation  by  TF-carrying  beads  (red 
fluorescence)  was  monitored  by  using  a  fluorogenic 
substrate  for  thrombin.  White  dotted  lines  indicate 
channel  walls  of  microfluidic  device.  (B)  A 
representative  TF-carrying  bead  did  not  initiate 

in  a  less  confining  channel  (160  x  160 
representative  TF-carrying  bead  did  initiate 
coagulation  (blue)  in  a  more  confining  channel  (10x10 
pm2).  Linescans  of  fluorescent  intensity  along  the  green 
dashed  line  are  given  below  the  images. 
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Figure  5.  In  2-dimentional  simulations  of 
sub-threshold  patches  of  TF  (1  nM; 
radius  r  =  3  pm)  placed  in  either 
confining  or  less  confining  channels  of 
plasma,  clotting  initiated  only  in  the 
confined  space.  (A)  A  schematic  drawing 
shows  the  device  design.  (B,C)  Plots 
show  [thrombin], otai,  the  sum  of  [thrombin] 
and  [meizothrombin],  including  both  free 
and  bound  species,  at  t  =  10  min.  (B)  In 
less  confining  (160  pm)  channel,  no 
coagulation  was  initiated  at  10  min,  and 
[thrombin], otai  remained  at  a  low  level 
(blue).  (C)  In  confining  space  (10  pm) 
containing  a  patch,  clotting  was  initiated 
and  [thrombin], 0,ai  increased  to  a  high 
concentration  (red).  No  clotting  was 
initiated  in  the  top  channel  (blue),  which 
did  not  contain  a  patch  of  TF  Scale  bars 
=  50  pm.  The  red  dots  in  the  schematic 
indicate  patches  of  TF,  and  they  are  not 
drawn  to  scale. 


Confinement  regulates  initiation  of  coagulation  by  bacteria  GMDs 

Having  tested  and  described  initiation  of  coagulation  by  TF  in  confined  volumes,  we  then 
tested  the  hypothesis  that  confinement  of  bacteria  could  lead  to  coagulation  of  plasma  in  vitro. 
We  tested  whether  a  sub-threshold  cluster  of  Bacillus  cereus  ( B .  cereus)  could  also  initiate 
coagulation  by  limiting  mass  transfer  of  active  factors  using  confinement.  In  these  experiments, 
we  placed  human  platelet  poor  plasma  with  sub-threshold  size  gel  microdroplets  (GMDs) 
carrying  B.  cereus ,  “bacteria  GMDs”,  in  large  microfluidic  chambers  with  diameter  of  9  mm  and 
height  of  1mm  and  in  small  microfluidic  channels  with  dimensions  40  x  40  pm*  (Fig.  6).  In  less 
confining  microfluidic  chambers,  active  factors  generated  by  individual  bacteria  GMDs  were 
constantly  removed  by  diffusion,  and  bacteria  GMDs  did  not  initiate  coagulation  for  over  20  min 
(Fig.  6 A):  however,  placing  similarly-sized  individual  bacteria  GMDs  in  confining  microfluidic 
channels  decreased  the  clot  time  to  less  than  10  min  (Fig.  6 B)  (N  =  3  beads  in  3  separate 
experiments,  p  <  0.05).  This  result  indicates  that  confining  microfluidic  channels  limit  the  mass 
transfer  of  active  factors  generated  by  bacteria  GMDs  and  allow  active  factors  to  accumulate 
locally  above  threshold  concentration  to  initiate  blood  coagulation. 
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Figure  6.  Bacteria  gel  microdroplets 
(GMDs)  initiate  clotting  rapidly  only 
when  confined.  Initiation  of 
coagulation  by  bacteria  GMDs 
(green  fluorescence)  was  monitored 
by  using  a  fluorogenic  substrate  for 
thrombin.  White  dotted  lines  indicate 
channel  walls  of  microfluidic  device. 
(A)  A  cluster  of  B.  cereus  on  a  GMD 
did  not  initiate  clotting  in  a  less 
confining  microfluidic  chamber  (9 
mm  diameter  x  1  mm  height)  for  20 
min.  (B)  A  cluster  of  B.  cereus  on  a 
GMD  did  initiate  clotting  (blue)  in  a 
confining  channel  in  10  min  (40  x  40 
pm2  channels). 


Chemical  Amplification  with  Threshold 

To  generate  a  threshold  we  decided  to  use  a  system  with  an  enzyme,  its  inhibitor,  and  its 
substrate,  in  which  the  input  is  the  enzyme  and  the  output  is  the  product  of  catalysis.  A 
theoretical  model  was  developed  to  predict  in  which  case  one  can  have  a  practical  threshold. 

The  existence  of  a  practical  threshold  is  generally  defined  as  a  big  relative  difference  between 
the  amount  of  output  when  the  input  is  slightly  below  the  threshold  and  that  when  the  input  is 
slightly  above  threshold.  Explicitly,  one  can  use  two  parameters  a  andy  to  define  the  existence 
of  practical  threshold  as 

( 1 )  Output  [ Input  =  (7  +a)  Threshold  J  -  Output  [  Input = ( 1  -  a )  Threshold ]  >  y  ( MaxOutput ) 


In  the  model,  by  letting  the  enzyme  equilibriate  with  the  inhibitor  and  subsequently  allowing  the 
unbound  enzyme  to  catalyze  the  production  of  the  product,  the  rate  of  production  of  the  output  is 


(2) 


So  +  KKf 

E available  =  —  ~  |  +  /Q  +  Kp  —  ^—\E0I()  +  (£'0+/0+A£))  j 


and  a  practical  threshold  exists  when 

y(ay  +  y-a) 

(3)  Threshold  =  I0>KD /V  f  r  1 

a  (1  -y) 


where  E0,I0,  and  S0  are  initial  concentration  of  the  enzyme,  inhibitor,  and  substrate. 


respectively,  kcal  and  Ku  are  Michaelis-Menten  rate  constants,  KD  is  the  dissociation  constant 
of  the  enzyme-inhibitor  complex. 


To  check  the  validity  of  the  above  results,  assays  were  performed  with  acetylcholinesterase 
(whose  inhibitor  is  m-isopropyl-2’,2\2’-trifluoroacetophenone,  with  KD  of  6.3  nM  as 

determined  by  simple  inhibitory  assays)  and  thrombin  (whose  inhibitor  is  hirudin  with  KD  of  20 

fM15).  Reaction  rates  were  normalized  by  the  rate  of  the  fastest  reaction  in  the  set.  Equation  (2) 
was  used  for  theoretical  data.  Experimental  initial  rates  were  calculated  from  time  trace  of 
fluorescence  intensity.  Indeed,  the  theoretical  model  and  experimental  data  agree  well  (Figure 

7). 

A  B 


Figure  7.  Agreement  between  theory  and  experiment  for  A)  acetylcholinesterase/  m- 
isopropyl-2’,2’,2’-trifluoroacetophenone  (1.7x1  O'7  M)  and  B)  thrombin/hirudin  (1.5x10‘9  M). 
Substrates  for  acetylcholinesterase  and  thrombin  were  indoxyl  acetate  (1.3x10‘sM)  and 
(Boc-Asp(OBzl)-Pro-Arg-MCA  (1.7x1  O'5  M),  respectively. 


Analog-to-digital  Conversion  of  Chemical  Signal  with  Application  for  Diagnostics 

Figure  8.  Threshold  with  AChE  (avidin-AChE 
conjugate)  with  visual  readout.  Inhibitor  concentration 
was  set  up  in  steps  of  1.3  fold.  A  1.3-fold  increase  in 
enzyme  concentration  is  detected  by  a  step  in  the 
position  of  the  threshold  (i.e.  an  increase  in  the  number 
of  completed  reactions).  The  reactions  took  place  in  a 
96-well  plate,  and  the  total  volume  in  each  well  was  43 
pL.  Concentrations  were  final  concentrations;  the 
concentration  of  ths  substrate  acetylthiocholine  was  0.93 
mM.  The  image  was  taken  at  39  minutes  with  a 
stereoscope  with  incident  light.  The  image  was 
subtracted  by  the  image  at  151  minutes  (when  all  wells 
had  reacted),  and  then  corrected  with  Photoshop 
(contrast+95). 

We  applied  chemical  amplification  with  threshold  to  converting  analog  chemical  signals  (e.g. 
concentration  of  a  substance)  into  digital  signals  (a  series  of  on/off  reactions)  that  describes  the 
analog  signal.  This  conversion  was  achieved  by  subjecting  the  input  concentration  to  different 
thresholds.  From  the  on/off  results  of  these  tests,  we  can  infer  quantitative  information  about  the 
input.  A  potential  use  of  this  technique  is  quantitation  in  resource-limited  settings  where 
traditional  techniques  such  as  absorption  or  fluorescence  require  unavailable  extra  equipment. 
Particularly,  we  applied  this  technique  into  an  immunoassay  for  cystatin  C,  a  biomarker  of 
kidney  function.  We  used  SlipChip16  as  the  microfluidic  platform  for  the  assay  because  it 


inhibitor 

3.44  4.47  5.81  7.56  nM 


completed  incomplete 


accommodates  multistep  processing  (required  for  the  immunoassay)  and  multiplexing  (required 
for  the  analog-digital  conversion). l7'19  The  chemistry  we  developed  uses  acetylcholinesterase 
(AChE)  as  the  reporter  for  the  immunoassay,  a  33-fM  inhibitor  to  set  up  the  threshold,20 
acetylthiocholine  as  the  substrate,  and  the  starch/I3"  complex  that  would  react  with  the  product  to 
become  clear  to  produce  a  visual  readout  (dark  to  clear).  We  first  tested  this  chemistry  in  well 
plates  (Figure  8).  This  chemistry  allows  for  visualization  of  the  signal  by  the  naked  eyes  and 
recording  the  results  with  a  consumer’s  camera.  We  were  able  to  distinguish  0.64  mg/L  cystatin 
C  (normal  concentration)  from  0.96  mg/L  cystatin  (stage  3  chronic  kidney  disease)21,2^  in  5  pL 
serum  as  shown  in  Figure  9.  The  average  difference  in  the  number  of  incomplete  reactions  was  3, 
and  the  95%-confidence  interval  was  0.9-5. 1. 

Figure  9.  Detection  of  1.5  fold  change  in 
concentration  of  cystatin  C.  (A)  Unmodified 
picture  of  SlipChip  (0.64  mg/L  trial  3  at  55 
minutes).  (B)  Distinguishing  between  0.64  mg/L 
cystatin  C  (normal  concentration)  and  0.96  mg/L 
cystatin  (stage  3  chronic  kidney  disease)21,22  in 
serum  diluted  2  times  in  PBS  buffer  (pH  7,  with 
Pluronic  FI  27  1  mg/mL).  5  pL  of  sample  in 
cystatin  C-spiked  serum  was  used  for  each  trial. 

Each  trial  required  3  one-hour  incubations: 
binding  with  the  antibodies  (one  of  which  bound 
to  beads  and  the  remaining  of  which  bound  to 
biotin)  followed  by  washing,  binding  with  avidin- 
acetylcholinesterase  followed  by  washing,  and 
reaction  with  the  inhibitor.  The  concentrations  of 


A 

4.5  mm 


the  inhibitor  droplets  are  shown  in  the  figure, 
while  the  concentration  of  acetylthiocholine  in 
the  substrate  mixture  droplet  was  1.0  mM.  The 
pictures  were  taken  after  55  minutes  of  the 
reaction  (except  for  0.96  mg/L  cystatin  C  trial  2 
which  is  presented  at  40  min,  at  which  point  the 
reactions  had  already  completed)  with  a  standard 
digital  SLR  camera  (Canon  RebelXT  with 
EFS18-55mm  lens).  The  images  were  rotated  and 
cropped  with  ImageJ,  then  adjusted  with 
Photoshop  (level  change  by  setting  a  spot 
between  wells  6  and  7  to  be  white, 
brightness +40,  contrast+ 1 00),  and  masked  to 
draw  attention  to  the  substrate  area  with 
Freehand. 
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Conclusions 

This  project  has  been  successful  in  demonstrating  that  we  can  begin  to  understand  the 
spatiotemporal  dynamics  of  complex  reaction  networks  by  using  a  chemical  model  of  a  network 
that  is  based  on  a  modular  mechanism. 

The  results  presented  here  demonstrate  how  local  confinement  of  a  small  amount  of 
stimulus,  such  as  TF  or  clot-inducing  bacteria,  can  dramatically  alter  the  outcome  of  coagulation, 
a  nonlinear  reaction  network  with  threshold  kinetics.  These  results  demonstrate  that  a  simple 
physical  phenomenon,  accumulation  of  activators  due  to  environmental  confinement,  can  have 
nontrivial  consequences  when  applied  to  a  reaction  system  having  threshold  kinetics.  We  also 
demonstrated  that  we  can  observe  a  chemical  threshold  with  an  enzyme  and  an  inhibitor.  In 


ongoing  work,  we  are  applying  this  threshold  to  detect  particles  (e.g.  bacteria)  in  microfluidic 
devices  with  stochastic  confinement.  Thresholds  abound  in  physiology,  including  thresholds  seen 
in  blood  coagulation2’-24  and  activation  thresholds  in  the  immune  system,2*  so  we  expect 
confinement  to  play  a  role  in  many  contexts.  In  ongoing  work,  we  are  investigating  the  use  of 
thresholds  for  chemical  amplification. 
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